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Within the Ginzburg-Landau approach a theoretical study is performed of 
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C3 , the effects of confinement on the transition to superconductivity for type-I 



and type-II materials with surface enhancement. The superconducting or- 
der parameter is characterized by a negative surface extrapolation length b. 
This leads to an increase of the critical field H c ^ and to a surface critical 



£>. \ temperature in zero field, T cs , which exceeds the bulk T c . When the sample 
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is mesoscopic of linear size L the surface induces superconductivity in the 
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superconductivity in thin films of type-I materials is akin to capillary conden- 



interior for T < T C (L), with T C (L) > T cs . In analogy with adsorbed fluids, 



sation and competes with the interface derealization or "wetting" transition. 
The finite-size scaling properties of capillary condensation in superconduc- 

c : 

Q , tors are scrutinized in the limit that the ratio of magnetic penetration depth 

O 

to superconducting coherence length, n = A/£, goes to zero, using analytic 

calculations. While standard finite-size scaling holds for the transition in non- 
zero magnetic field H, an anomalous critical-point shift is found for H = 0. 
The increase of T c for H = is calculated for mesoscopic films, cylindrical 
wires, and spherical grains of type-I and type-II materials. Surface curva- 
ture is shown to induce a significant increase of T c , characterized by a shift 
T C (R) — T c (oo) inversely proportional to the radius R. 



I. INTRODUCTION 

The Ginzburg-Landau (GL) theory of superconductivity continues to deliver surprises 
[1]. In this paper we focus on some remarkable consequences of special boundary conditions 
that enhance superconductivity at the surface of the material. In fact, superconductivity 
is already known to be enhanced for the common situation of surfaces against vacuum or 
insulators, as was demonstrated by the discovery of the surface critical field H c % [2]. We 
consider, however, different surfaces that enhance superconductivity more strongly. Within 
GL theory this is embodied phenomenologically by taking the surface extrapolation length 
b to be negative. It was shown that this not only leads to a further increase of the surface 
critical field H& [3], but also to an increase of the surface critical temperature in zero field, 
T cs [4]. The simple relation £,(T CS ) = —b, with £(T) the superconducting coherence length 
in bulk and in zero field, governs the shift from T c to T cs . Furthermore, for b < interface 
derealization transitions, which are the precise analogues of wetting transitions in adsorbed 
fluids, have been predicted for type-I superconductors [5]. 

For a semi-infinite system with a planar surface the GL surface free energy functional, 
including the boundary contribution, reads 

7[M = ^I#>)| 2 (1.1) 
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The magnetic field is taken parallel to the surface. For this orientation the interface dereal- 
ization or wetting transition can occur, provided k = A/£ < l/\/2 (type-I superconductors) 
and b < 0. Here, A is the magnetic penetration depth and £ is the zero-field superconducting 
coherence length. For k < 0.374 the wetting transition is of first order, and is accompanied 
by a prewetting line that extends into the bulk normal phase in the H — T phase diagram, 
and terminates in zero field at T cs [5]. An experimental realization of the prewetting phe- 
nomenon is, in hindsight, provided by the twinning-plane superconductivity transition in Sn 
[4]. The wetting transition itself has so far not been verified experimentally. For k > 0.374 



the wetting transition is predicted to be critical, without a prewetting line [5]. 

In this paper we study the effect of confinement on the wetting phase diagram, and, 
in particular, we examine the increase of T c . The situation we consider is analogous to 
that of a fluid adsorbed between parallel walls, which undergoes capillary condensation [6]. 
This phenomenon occurs slightly below the saturated vapour pressure, and arises from a 
competition between surface contributions to the free energy, proportional to the surface 
area, and volume contributions, proportional to area times wall separation L. For large 
L, the pressure or chemical potential for which the fluid condenses between the walls, is 
shifted by a small amount, proportional to 1/L, from the usual bulk coexistence line. In the 
presence of a wetting transition for the semi-infinite system, there is an interesting interplay 
between capillary condensation and prewetting, leading to surface triple points. We study 
the counterparts of these phenomena for type-I superconductors, in the \ow-K limit. 

In zero field, the increase of the surface critical temperature T cs for samples with b < is 
not limited to type-I materials, but occurs for type-I and type-II alike. It is therefore justified 
to devote special attention to the effect of confinement on this phenomenon. The increase 
of T c is unique to superconductivity, since in fluids confinement generally suppresses the 
critical point of phase separation. In contrast, we find that in superconductors the critical 
temperature of a mesoscopic sample with surface enhancement not only exceeds the bulk 
T c , but is also greater than T cs . We study this effect for planar films, cylindrical wires and 
spherical grains. For surfaces with curvature an important additional increase of T c is found. 

The assumption b < is crucial here, since for b = oo (surfaces against vacuum or 
insulators) and a fortiori for b > (surfaces against normal metals or ferromagnets) there 
is no increase of T c relative to the bulk value. For b = oo the effect of confinement leads to 
a well-documented increase of the critical field and the presence of a tricritical point where 
the transition to superconductivity changes from second-order to first-order as the field is 
increased. The effects of sample topology for this case (b = oo) have been the subject of 
thorough experimental [7] and theoretical [8] investigation. 

This paper is organized as follows. In Section 2 we study the effect of confinement 



on the wetting phase diagram. The limit of strongly type-I superconductors turns out to 
be very instructive here, since various analytic results can be obtained for k — > 0, the 
details of which are outlined in Section 3. The finite-size scaling properties of the capillary 
condensation transition in non-zero field, and the link to the anomalous critical-point shift 
in zero field are addressed here. In Section 4 we derive and discuss the increase of T c for 
mesoscopic surface-enhanced superconductors. Conclusions and remarks pertaining to the 
experimental relevance of our results are presented in Section 5. 

II. CAPILLARY CONDENSATION AND PREWETTING FOR STRONGLY 

TYPE-I SUPERCONDUCTORS 

In this Section we discuss the precise analogy between the capillary condensation tran- 
sition in a fluid confined between parallel walls and the transition to superconductivity of a 
mesoscopic film of type-I material in a parallel magnetic field. The surfaces of the film are 
characterized by surface enhancement of superconductivity (negative extrapolation length 
b) and we consider the case of identical surfaces, which is sufficient to address the basic 
phenomena. For temperatures sufficiently close to T c , interface derealization comes into 
play and allows us to study, in close analogy to what may happen in a confined fluid, how 
capillary condensation competes with the prewetting phenomenon. 

There are four relevant lengths in our system: the magnetic penetration depth A, the 
coherence length £, the surface extrapolation length b, and the film thickness L. In order 
to study the interplay between capillary condensation and prewetting most clearly and 
accurately, it is very useful to take the limit k = A/£ — > 0, corresponding to extreme type-I 
superconductors. It is important to specify that in taking this limit, we let A tend to zero, 
while keeping the other three lengths finite. In this limit not only is the wetting transition 
of first order but also the prewetting transition remains of first order down to zero magnetic 
field, so that the competition with the capillary condensation transition (also of first order) is 
not complicated by second-order nucleation phenomena that occur for k > 0. Furthermore, 



the vortex phase, which we find to play a role even for k considerably less than l/\/2 in a 
film with enhanced surfaces, is unimportant at k — 0. 

Besides these reasons pertaining to clarity, the limit k — ► offers the major advantage 
that the problem can be studied analytically, and the important finite-size scaling laws for 
the phase transitions can be calculated exactly (see next Section for details). Several of these 
laws continue to hold, in as far as the leading singularity in the asymptotic regime of thick 
films is concerned, for small k > 0, as long as the phase diagram undergoes quantitative 
changes only. Therefore, many of the properties that we can demonstrate analytically at 
K — 0, serve as a good first approximation for a significant part of the type-I regime. We 
have verified this by numerical computations for k > 0. 

The usefulness of taking the zero-K limit has already become clear in previous studies of 
interfacial properties in type-I superconductors, most notably in the derivation of an exact 
interface potential for wetting and prewetting [9]. Moreover, it has been shown extensively 
that the thermodynamic behaviour at k > (but not exceeding l/\/2) can often be captured 
by means of rapidly converging expansions in the parameter k [10,11]. 

The two basic physical states of the film consist of either superconducting surface sheaths, 
extending from one or both surfaces into the interior, or a superconducting film state, which 
occupies the whole space between the surfaces. The former correspond to prewetting layers 
and the latter to capillary condensation. For computing these states we recall that for 
k = the magnetic induction A(x) and the superconducting wave function ip(x) exclude 
one another in space [9]. Furthermore, since A is a simple step function, the only pertinent 
GL equation is that for if), which after suitable rescaling (as in [9]) reads 

^ = ±^ + ^ 3 (II.l) 

The +(— ) signs pertain to T > T C {T < T c ), respectively. Note that, with the present 
rescaling convention, lengths are measured in units of the (zero-field) coherence length £, 
and in the bulk superconducting phase, ip — 1, 
The boundary conditions are 



m = ^(0)/6 (II.2) 

The useful first integral of (II. 1), 

^ 2 = ±ip 2 + ^/2 + C (II.3) 

allows one to employ a simple phase-portrait analysis for determining the characteristics 
of all possible solutions. The integration constant C is determined using the boundary 
conditions. For T < T c capillary condensation states exist for C < 1/2, while for T > T c 
they occur for C < 0. They are symmetric with respect to the middle plane of the film, 
ijj(x) = ip{L/^ — x). For these states ip(x) has the shape of a "hammock", with a minimum 
at x — L/2£. C is a smooth function of L which tends to 1/2 for large L, as ip{L/2^) tends 
to the bulk value 1. An interesting point to note is that the magnetic field H is fully expelled 
in these states and therefore the profiles ip(x) do not depend on H . In particular, C(L) is 
independent of H. 

In contrast, the prewetting states depend on the applied field. These states are charac- 
terized by profiles ip{x), interrupted by a magnetic "gap" in which ip(x) = 0. The phase 
portrait analysis indicates that two types of solutions can be considered: symmetric states 
consisting of two superconducting surface sheaths located oni6 [0, //£] and \{L — /)/£, L/£\, 
with a central gap separating them, or asymmetric states with a sheath at one surface only, 
on x G [0,//£], followed by a gap extending to the other surface. In practice the asym- 
metric states are irrelevant, even for I > L/2, since their free energy is higher than that 
of symmetric states (either of prewetting type, with a central gap, or of capillary conden- 
sation type, without a gap). For prewetting states, C depends on H and not on L. Its 
form, C(Hr) = Hft, is the same as for the semi-infinite system, since the "gap" acts in the 
same way as a normal phase (N) in bulk. The quantity Hr is a reduced field defined in [9] 
and related to H in the manner Hr oc £?H . The value of H\ determines the magnitude 
(squared) of the gradient of ip at the interior points x = //£ and x — (L — /)/£ at which ip 
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Ptp C ap(0) 

L = 2U #(±^ 2 + ^ 4 /2 + C)- 1/2 (II.4) 

Jlbm. 



vanishes. For T < T c , at bulk two-phase (SC/N) coexistence, C = 1/2, while in the bulk 
normal phase, C > 1/2. The term "bulk" refers to an infinite system (without surfaces, in 
principle). 

We can conveniently express the film thickness using the profile of a capillary condensed 
state ("cap"), through the relation 

Ari>ca P (0) 

Here, ip m = ?/>(L/2£), the value in the middle of the film. Likewise, we can obtain the 
thickness of a surface sheath in the prewetting state ( "pw" ) through 

rippw(0) 

l = U d^j(±^ 2 + ij 4 /2 + Hi)- 1 ' 2 (II.5) 

Jo 

For explicit expressions for ip m , ipca P (0) and ip pw (0), see Appendix A. 

Similar compact expressions are available for the reduced (i.e., dimensionless) free ener- 
gies. For capillary condensed states, 

ArV'cap(O) 

Using (II.4) this can be simplified by separating out the dependence on the magnetic field, 
which is just H^L/l;, the free energy cost of expelling the field over the whole thickness of 
the film. 

For prewetting states the free energy of a symmetric state with two surface sheaths is 

rippw(0) 

lpw = 2 dmi - rm(±r + r/2 + #ir i/2 (n.7) 

Jo 

In order to be able to discuss the phase diagram for temperatures below, above, and at 
bulk T c , it is convenient to express the thickness of the film in units of |6| instead of £ (since 
£ diverges at T c in zero field). It is understood that the value of b is the result of the surface 
preparation of the sample, and can therefore be considered a material constant within the 
explored ranges of field and temperature. 

In order to show most clearly the topology of the new phase diagram of the thin film 
with surface enhancement, we have chosen the (reduced) thickness L/\b\ = 8. The result 
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leap = 2 / #(tf| - ^ 4 /2)(±^ 2 + V>V2 + CT 1/2 (H.6) 



is presented in Figure 1. The temperature variable is t = (T — T C )/(T CS — T c ), so that the 
first-order interface derealization transition, or "wetting" transition, is located at to < 0, 
while the bulk critical point in zero field is at t c — and the surface critical point in zero 
field is at t cs — 1. The magnetic field H is in units of Hp, the wetting transition field. The 
ratio H/H D is related to H R through the equation H/H D = \/2H R {S,/b) 2 D /{^/b) 2 . The thin 
straight line from D to the origin is the bulk two-phase coexistence line. The new phase 
transitions relevant to the mesoscopic film are indicated by thick solid lines. 

The main transition is the capillary condensation line, which consists of three parts. 
For high fields this line is more or less parallel to the bulk reference line, and represents a 
transition from a fully normal film to a fully superconducting film. Between Tl and T2, 
however, for decreasing H, capillary condensation is preceded by the prewetting transition. 
The film thus goes superconducting in two distinct steps: (i) from a fully normal state to a 
state with two superconducting surface sheaths and a normal gap, and (ii) from the latter 
to a fully superconducting film. At transition (ii) the gap between the surface sheaths is 
still finite. Incidentally, we can compute the line in the phase diagram on which I = L/2, so 
that the gap vanishes and the two surface sheaths touch one another. For all temperatures 
between the wetting point and the prewetting critical point, this line lies at lower fields 
than the capillary condensation transition, and consequently has no physical significance. 
Finally, for temperatures between that of T2 and t c (L), the transition proceeds in a single 
step, from fully normal to fully superconducting. 

The prewetting phenomenon is thus confined to an "island" in the phase diagram, where 
a film with superconducting surface layers and a normal interior is thermodynamically stable. 
All transitions in non-zero field are of first-order. The points Tl and T2 are genuine triple 
points of the film. The three coexisting film phases are represented by their wave function 
profiles in Figure 2. Likewise, Figure 3 illustrates the triple point T2. The prewetting line 
(between Tl and T2) lies exactly on the prewetting line of the semi-infinite system (dashed 
line), which extends from the wetting point D to the surface critical point of the semi-infinite 
system in zero field, at t — t cs — 1- 
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In zero field the capillary condensation ends in a critical point, at t — t c (L). This critical 
point will be discussed in detail in Section 5. We shall derive there that t c (L) is only slightly 
above t cs . On the scale of the figure the two points appear coincident. 

Upon lowering L/\b\ the points Tl and T2 approach each other, and for L/\b\ between 
6 and 7 the prewetting "island" vanishes. For L/\b\ < 6 only the capillary condensation 
transition remains. For example, for L/\b\ = 1 the capillary condensation line appears as 
a straight line parallel to the bulk coexistence line, and ends in zero field at t c ~ 2.4 (cf. 
Section 4). 

On the other hand, upon increasing L/\b\ the triple point T2 moves rapidly to zero field 
and t — t cs , an d Tl moves slowly towards the wetting point D. The prewetting line remains 
fixed. The capillary condensation transition converges, for T < T c , to the bulk coexistence 
line. However, for T > T c , the capillary condensation line converges to the segment [T c , T cs ] 
of the temperature axis (at H = 0). This is a consequence of the anomalous critical-point 
shift in zero field. The phase boundary thus develops a corner singularity at the origin 
(t — 0, H — 0). The precise manner in which the phase boundary scales in the limit L — > oo 
is the subject of the next Section. 

III. FINITE-SIZE SCALING OF CAPILLARY CONDENSATION 

In order to examine how the capillary condensation phase boundary approaches the bulk 
coexistence line in the limit L — ► oo we distinguish the following regimes. 

3.1. T < Tjj: below the wetting transition, approaching partial wetting. 

In this regime the complication of surface superconductivity does not arise and the transi- 
tion is from the normal phase directly to a superconducting film with complete expulsion of 
the magnetic field. The transition occurs when the free energy 7 cap , given in (II. 6), equals 
that of the normal phase, which is zero. For large L this condition is very well approxi- 
mated by replacing the upper and lower limits of the integral by their asymptotic values, 
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V'cap(O) — > i/)(0) an d V'm - ► 1- Here, ^(0) is the surface value of the wave function profile 
associated with the superconducting phase in bulk, at temperature T. This leads to the 
familiar result, akin to Laplace's or Kelvin's equation for a confined fluid [6], expressing the 
free energy balance between a cost in bulk and a cost in surface contributions, 

{H 2 R - l/2)L/£ = -2 lw , sc (III.l) 

The r.h.s. is by definition — lim£_ >0O 7 cap and represents (minus) the surface free energy 
of two wall/SC "interfaces". Since ^w,n = 0, in the absence of superconducting surface 
sheaths, the r.h.s. actually equals 2'jsc/N cos 9, familiar in the context of Young's equation 
for the contact angle 9 in the partial wetting regime. The l.h.s. gives the net free energy 
cost, per unit volume, of expelling the magnetic field (cost H R ) and, simultaneously, going 
superconducting (gain 1/2), multiplied by the thickness of the film. This net cost is positive 
for fields higher than the coexistence field (commonly referred to as critical field of the 
superconductor) H Rc = \j\pl. Equation (III.l) predicts that the capillary transition field 
Hji(L) approaches the coexistence field H Rc according to the power law 

(H R (L) - H R , C ) ~ v^^f^ (III.2) 

The exponent of L, — 1, simply reflects the difference between the surface dimension d — 1 
and the bulk dimension d. 

Numerical computations show that (III.l) is extremely accurate, even for thin films. For 
example, for "temperature" £/6 = —0.5, (III.l) is satisfied to an accuracy of 0.1% already 
for L/\b\ = 2. The deviation for L/\b\ = 1 is about 3%. Therefore, for practical purposes, 
(III.l) is correct for L/\b\ > 2. 

The wetting transition (D) occurs at £/6 = (£/6)d ~ —0.60, and for temperatures 
Td < T < T c (III.l) must be modified as follows. 

3.2. Td < T < T c : the prewetting regime, approaching complete wetting. 

In this regime the capillary condensation competes with the prewetting transition. To a first 
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approximation the surface free energy balance takes the Laplace or Kelvin form, analogous 
to (III.l), 

{H 2 R - l/2)L/i = 2( lw , N - lw , sc ) = 2 lsc . N (111.3) 

where ^w,n is the surface free energy of a semi-infinite system with a macroscopic surface 
superconducting layer (complete wetting). The last equality expresses that Antonov's rule 
holds for complete wetting [12]. 

However, this approximation is too crude. It neglects first of all that L should be replaced 
by L — 21 to take into account the thickness of the prewetting layers, constituting the part 
of the film which is already superconducting before capillary condensation occurs. But even 
with this correction, the resulting approximation is still not satisfactory, in comparison with 
numerical computations. In what follows we derive an accurate analytic approximation for 
large L. 

We start from the exact condition for capillary condensation 

ipw leap V / 

The magnetic field terms in these free energies lead to contributions 2H R l/£ and H R L/^, 
respectively, as is seen from the set of equations (II.4)-(IL7). In the limit L — > oo, the 
capillary condensation field H R approaches the bulk coexistence field H Rc = l/v^2, and the 
prewetting layer thickness / diverges, but very slowly. To see this in detail, we work out the 
integral and obtain the magnetic-field dependence of I, 

l{H R )/i = -^ In ff2 \ + k/i + o(l) (III.5) 

The divergence is only logarithmic, so that the constant /i/£ is an important correction for 
numerical purposes. Furthermore, it is the only important correction, since we verified that 
the remainder is insignificant, up till H R — 1/2 = 0(1). We remark that the upper spinodal 
of the prewetting transition occurs at H R = (1 + (£/fe) 2 ) 2 /2. The remainder is indicated by 
o(l), which signifies that it goes to zero as H R — > 1/2. Numerically these terms are found 
to vanish as H R - 1/2, or (H R - 1/2) ln(l/(ifjj, - 1/2)). 
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The constant h/C, can be calculated analytically, with the result 

/i/e = v / 21n2 + ^ln ^S°l~ 1 +2 f°° du{2u 2 + ir 1/2 - 2 H du{2u 2 )' 1/2 (III.6) 
V2 V(0) + 1 "'o J\ 

The value of ^(0) here corresponds to the limit of bulk two-phase coexistence, and is deter- 
mined through t/>(0) 2 = 1 + (i/b) 2 + ((1 + {£/b) 2 ) 2 — l) 1 / 2 . Typical values of Zi/£ are of order 
1, confirming the importance of this constant next to the leading logarithm in (III. 5). For 
example, for £/6 = —1, li/£ m 1.640. 

The geometrical interpretation of /i/£ is straightforward. Keeping only the leading and 
next-to-leading terms in 1{H^)/C, we arrive at the identification, 

h/i » KH^/Z, (III.7) 

with H\ x — \j1 = 1. This is qualitatively correct. For example, for "temperature" £/6 = — 1, 
li/£ = 1.640 while l{Hn t i)/£ = 1.504. Thus, li/£ corresponds essentially to the thickness 
of a thin surface sheath at a magnetic field well above the critical field. This thickness (in 
units of £) is of order 1. Consequently, the leading logarithm in (III. 5) gives the intrinsic or 
"net" thickness of the wetting layer which develops close to bulk coexistence. 

Having established the slow divergence of /, and contrasting it with the more rapid 
divergence of L, which is essentially proportional to 1/{H\ — 1/2), we collect carefully all 
terms proportional to H\ — 1/2 in (III.4) and find 

{H 2 R - 1/2)(L - 21) H = 2 1s c,n + {H 2 R - 1/2)(V2 + o(l)) + 0{e~^), (III.8) 

where o(l) vanishes for Hi — ► 1/2. A summary of the derivation is given in Appendix 
A. The l.h.s. features the net cost of expelling the magnetic field, while the first term 
on the r.h.s. gives the cost of having two SC/N interfaces. Especially interesting, and 
calculable analytically, is the correction to the surface tension, {Hi — 1/2) \/2, appearing as 
the second term in the r.h.s. Precisely in view of the slow divergence of / this contribution is 
numerically significant in combination with the l.h.s. Taking it into account greatly improves 
the accuracy of the approximation. 
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The surface tension correction has an interesting physical interpretation. In the complete 
wetting regime at bulk two-phase coexistence a superconducting/normal (SC/N) interface 
constrained at a distance / from the surface has a free energy (per unit area) that is higher 
than that of an equilibrium interface (infinitely) far away from the surface, by an amount 
which is given by the so-called interface potential V(l). This excess free energy is known 
exactly in the k = limit [9], and we are concerned here with the tail of V(l) for large /, given 
by V(l) oc exp(— V2V0- Therefore, the free energy cost of a constrained interface is easily 
found, by inserting the logarithmic divergence (III. 5), to be proportional to H\ — 1/2. Thus 
we arrive at the interpretation that the surface tension correction is due to the interaction 
or "interference" of the interface with the surface, which in the broader context of confined 
interfaces is often referred to as entropic repulsion [13]. 

In the same spirit, a correction is present for the capillary condensed superconducting 
state, relative to a superconducting state with infinite surface separation L. This correction is 
due to the interaction between the surfaces bounding the film, and also decays exponentially 
with separation. However, since the distance is now L instead of /, this contribution is of 
order exp(— v2£/0> which is negligible for our purposes (see Appendix A for details). 

In conclusion, in the complete wetting regime the finite-size shift of the transition to 
superconductivity has the same asymptotic scaling behaviour Hr — H^ c oc 1/L as in the 
partial wetting regime, but quantitatively L must be shortened by twice the wetting layer 
thickness I, and an effective further correction a/2£ must be subtracted from L in order to 
take into account the distortion of the two constrained SC/N interfaces. 

3.3. T = T c : the bulk critical isotherm. 

The finite-size scaling properties of the transition to film superconductivity at T c are in- 
teresting and merit a separate study, since they invoke universal quantities associated with 
the bulk critical point. At T = T c the zero- field coherence length is infinite, and we cannot 
use it as the unit of length. Instead we use |6|. The wave function must also be rescaled, 
because the normalization ipbuik = 1 is inconvenient at T c . Simple universal GL equations 
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result when we rescale x — ► (£/\b\)x, ip — > (^/|&|)V' — 0; -Hfl ~~ * (£,/b) 2 H R = Hr. The ratio 
H/Hd is invariant and equals \f2h,R(£,/b) 2 D . 
The GL equation now reads 



(III.9) 



and the boundary conditions take the form 



0(0) = -0(0) (111.10) 

0(L/|6|) = 0(L/|6|) 

Writing the first integral of (III. 9) as 

2 = 4 /2 + c (III.ll) 

we obtain for the film thickness, 

L = 2\b\ d<j){<j) 4 /2 + cr 1/2 , (111.12) 



with 4> m = (—2c) 1 / 4 and mp (O) 2 = 1 + (1 — 2c) 1 / 2 . For L — ► oo, c approaches zero from 
below and <p m vanishes. For the thickness of the superconducting surface sheath we have 

l=\b\ #(0 4 /2 + 4)" 1 / 2 , (111.13) 

Jo 

with <f) pw (o) 2 = i + (i - 2h R y/ 2 . 

We remark that, although T = T c marks the terminus of bulk two-phase coexistence, 
two-phase coexistence for the mesoscopic film continues to exist. Therefore, we will continue 
to use the terminology "capillary condensation" and "prewetting" in the same sense as in 
the previous subsections. 

Before discussing the free energies we examine how I behaves when the field h is turned 
to zero. A simple rescaling in (III. 13) suffices to extract the leading term, 

l/\b\ w h R 1/2 I" dx(l + x A /2)- l/2 w 2.20488 h R 1/2 (111.14) 

Jo 
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This power-law divergence is much faster than the logarithmic behaviour found in the prewet- 
ting regime below T c , approaching bulk two-phase coexistence. Experimentally, this implies 
that the diamagnetic response due to the surface superconducting sheath may be easier to 
detect when lowering H at T = T c than at T < T c . 

A similar reasoning leads to a simple relation between L and c, in the thick film limit, 

L/\b\ « (-c)- 1 / 4 2 5 / 4 / dx{x A - 1)- 1/2 (111.15) 

The integral equals 1.31103. So we conclude that c(L) decays as a power law, in contrast 
with the exponential decay seen in Appendix A, Eq.12. 

We now turn to the free energies. For capillary condensed states at bulk T c , 

lca P = 2 / rf0(4-0 4 /2)(0 4 /2 + c)- 1 / 2 , (111.16) 

J<t>m 

while for prewetting states, 

lpw = 2 / dd>(h 2 R - 4 /2)(0 4 /2 + h 2 R y l '\ (111.17) 

Working out the condition j pw = 7 cap for capillary condensation we find, 

h\{L - 2l)/\b\ = 5 lpw - 5 lcap (111.18) 

Using (III. 14) we see that the second term in the l.h.s. is of order h R . The first term on the 
r.h.s. is the surface free energy cost of constraining a surface sheath at H = and T = T c to 
terminate at x — l/\b\ instead of assuming its equilibrium power-law decay <fi{x) oc 1/x. This 
power-law decay is the analogue of "critical adsorption" for fluids [14]. Analytic calculation 
gives 

hpw = h¥V2 f°° dxx 2 {l - (1 + 2/rr 4 )- 1 / 2 ) (111.19) 

The integral equals 1.03939. The constrained surface sheath can be interpreted as a con- 
strained interface interacting with the surface. This interpretation is quite unconventional 
in this case, since an equilibrium interface does not exist at bulk T c . Nevertheless, assuming 
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the existence of an interface potential V(l) for the constrained interface leads us to infer 
V(l) oc l~ 3 , in view of (III. 14) and (III. 19). The exponent —3 is reminiscent of finite-size 
interactions at bulk criticality and will show up again in the next paragraph. Incidentally, 
note that 5^ pw > 0, corresponding to repulsion between the surface and the constrained 
interface (cf. "unlike" surfaces repel). 

The second correction <57 cap is the finite-size interference free energy between the two 
surfaces at separation L bounding the film, in the superconducting state at bulk T c and 
in zero field. This interference is akin to the generalized Casimir effect [15,16]. Analytic 
calculation gives 

hca P = -2 5/4 (-C) 3 / 4 (^^^(l - 1/0T 1 / 2 - jH<W) (HI.20) 

The integrals add up to 0.43701. Converting the |c|-dependence into an L-dependence, using 
(III. 15) we observe that the finite-size interaction is attractive ("like" surfaces attract) and 
decays in the manner L~ 3 . We scrutinize this generalized Casimir effect for superconductors 
and the associated universal exponents and amplitudes elsewhere [17]. 

Returning now to (III. 18) we see by simple inspection that h,R oc L~ 2 , and that all the 
leading corrections we extracted are of the same order, L~ 3 . We are therefore left with the 
simple problem of solving for the amplitude A in the asymptotic behaviour 

h R *A(L/\b\)- 2 (111.21) 

Numerical solution gives A rj 36.2869. The fact that the exponent of L equals —2 is linked to 
the fact that the mean-field value of the critical exponent vu is 1/2. This exponent describes 
the divergence of the field-dependent coherence length £(/ir) along the bulk critical isotherm 
approaching the bulk critical point, 

i(h R ) oc h R v » (111.22) 

The transition to superconductivity for the film occurs when £(/ir) ~ L, whence (III. 21). 
The sense in which universality holds is governed here by the validity of mean-field theory 
for classical superconductors. 
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A numerical computation of the finite-size shift of the critical field at bulk T c supports 
the analytic leading result (III. 21) and suggests that the next-to-leading term is of order 
(L/\b\)~ 3 , implying slow convergence. For L/\b\ = 10 the correction to the leading term is 
about 31%, while at L/\b\ = 100 the correction is about 3.7%. 

3.4. T c <T < T C (L): the bulk supercritical region. 

Even though T > T c , in this regime we still find a competition between prewetting-like states 
and capillary condensation. The main difference with respect to the prewetting region below 
T c is that an equilibrium SC/N interface and hence also its surface tension no longer exist. 
Therefore, the main modification to (III. 8) is that the first term on the r.h.s. is absent for 
T > T c . Furthermore, since the critical field is zero, the combination H\ — 1/2 simplifies to 
H\. The other modifications to (III. 8) will now be studied in detail. 

We start, as usual, from (III. 4). The magnetic field terms again lead to the net free energy 
contribution H\{L / '£ — 2/ / '£) . The prewetting layer thickness I has a different interpretation 
than for T < T c . Above T c no thick wetting layer can develop, since the infinite system 
consists of a single normal phase only. So, I just measures the extent of penetration into 
the bulk of the tail of the surface superconducting sheath. As the field H goes to zero the 
superconducting wave function no longer vanishes at x — l/C, but decays exponentially as 
a function of the distance x from the surface, so that, mathematically, / diverges although 
physically the penetration is of short range only. 

To see how I behaves above T c , we employ (II. 5) with the + sign. We obtain 

l(H R )/S = In -i- + hH + o(l) (111.23) 

tlR 

This is similar to (III. 5). The constant ^/C can be calculated analytically, with the result 

h/( = hl #(£)_ (m24) 

1 + \/l + </>(0) 2 /2 
The value of ^(0) here corresponds to the zero-field limit, and is determined through ^(O) 2 = 
2(— 1 + (£/6) 2 ). Typical values of /2/C are of order 1. For example, for £/6 = — 2, ^/^ « 1.184. 
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However, Z 2 /£ approaches zero and changes sign as £/|6| is decreased to about 1.07, not far 
from the surface critical point T cs . 

The geometrical interpretation of / 2 /£ is similar to that of Zi/£ discussed previously for 
T < T c . To a first approximation we can identify 

k/S » /(Hfl,2)/^, (HI.25) 

with H R2 — 1- This is reasonable. For instance, for £/6 = —2, Z 2 /£ = 1-184 while 
Z(i?# ;2 )/£ = 1.381. So we arrive at the interpretation that £ 2 /£ corresponds to the thickness 
of a t/im surface sheath in a finite field (of reduced strength unity). This interpretation 
can only be used as long as the (reduced) spinodal field exceeds 1. The spinodal line for 
prewetting states above T c is determined by H R = ((£,/b) 2 — \)j\[2. Consequently, (III. 25) 
makes sense as long as £,/\b\ > 1.554. The remainder o(l) in (III. 23) appears to vanish in 
the manner H R \n(l/H R ) as follows from numerical inspection. 

We now return to the condition for capillary condensation, which can be written as 

H&L - 21)/£ = r^ dM\tf + ^ 4 /2 + C)- 1 / 2 - r M0) dM\l? + ^ 4 /2 + H 2 R )- 1 ' 2 

(111.26) 

with i, 2 m = -1 + (1 - 2C) 1 /2, ^ cap( o)2 = -1 + ^/bf + ((-1 + {i/b) 2 ) 2 - 2C) 1 / 2 , and 
iP pw (0) 2 = -1 + (^/6) 2 + ((-1 + (i/b) 2 ) 2 - 2H 2 R ) l l 2 . The two integrals can be studied fairly 
easily, since expanding in the small parameters H R or \C\ poses no problems regarding 
the exchange of differentiation and integration, in contrast with the case T < T c . We 
find analytically that the integrals are, in leading order, simply proportional to C and H R , 
respectively. The result is 

H 2 R (L - 21) /I- = H 2 R (l + o(l)) - C(l + o(l)) (111.27) 

where the two terms o(l) vanish in the limits H R — > and C — > 0, respectively. 

We are thus left with the final task of determining the dependence C(L). This can also 
be done analytically, starting from (II. 4). In the limit C — * 0, with C < 0, we readily find 
the leading behaviour, 



LIZ = H-l/C) + L 2 /Z + o(l) (111.28) 

The constant is given by the expression 

j = 2/ 2 /e, (HI.29) 

with I2/Z as given in (III. 24). This quantity varies smoothly between —00 for £/6 = — 1 
(surface critical point T cs ) and the value 5 In 2 m 3.466 for Z/b — > —00 (bulk critical point 
T c ). It changes sign at "temperature" £/6 ~ —1.064. The remainder o(l) is numerically 
found to be proportional to Cln(— 1/C). 

Inverting (III. 28) to get C(L) we arrive at the following conclusion, which is the coun- 
terpart of (III. 8) for temperatures above T c , 

H 2 R (L - 21) /Z = H 2 R (1 + o(l)) + e L ^ e~ L ^ (111.30) 

which implies an exponentially fast decrease of the field as a function of L or, equivalently, a 
logarithmic divergence of L as a function of 1/Hr. This is in sharp contrast with the simple 
power-law found for the usual capillary condensation below T c . 

Since now I and L behave essentially in the same manner (logarithmic) as a function of 
1/Hr, we investigate numerically the interesting difference L — 21, in the limit H — > 0. Using 
(III. 23), (III. 28) and taking advantage of the equality (III. 29), we find that the constants 
cancel and we are left with 

(L - 2l)/i « H-Hl/C) (111.31) 

On the other hand, (III. 30) implies 

LIZ - 21 IZ - 1 » -C/Hl (111.32) 

In combination with the previous result the difference A = {L—2l)/Z must solve the equation 

A = ln x ^ T (111.33) 

Numerically, this gives A ps 1.2785. In conclusion, the difference L — 21 converges to a 
finite length, as we follow the capillary condensation transition into the asymptotic regime 
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L — > oo. We have verified this analytic result numerically, and the agreement is very good 
for sufficiently large L. For instance, for C,/b = —1.5, A is reproduced to 4 digits if we take 
L/\b\ = 20. For £/6 = —2 we achieve similar accuracy taking L/\b\ = 25. 

IV. CRITICAL-POINT SHIFT IN ZERO FIELD 

In zero magnetic field the transition to superconductivity is of second order and can be 
calculated using the linearized GL equation. This can be seen by calculating, in the pres- 
ence of a magnetic field, the location of the tricritical point where the order of the transition 
changes from second to first order, as the field is increased. Furthermore, in zero field the 
dependence on the GL parameter k drops out so that the results are valid for all classical 
superconductors, regardless of their type. We study three different geometries: planar slab, 
cylindrical wire and spherical grain. For each geometry we calculate the critical temperature 
as a function of the thickness (or diameter) of the mesoscopic system with surface enhance- 
ment of superconductivity. 

4.1. Planar film. 

For the planar sample, a slab or film with two parallel surfaces, we allow in general a dif- 
ferent enhancement on each surface. Thus we assume two surface extrapolation lengths, b\ 
and &2- Scaling all lengths with the zero-field coherence length £ leads to the GL equation 

i> = ±ip (iv. i) 

with boundary conditions 

^(0) = (^/6i)V(0) 

^(L/0 = -(£/6 2 )^(L/0 (IV.2) 

Since we are interested mostly in enhancing superconductivity (b < 0) we are concerned 
with T > T c , corresponding to the + sign in the r.h.s. of (IV. 1). 
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Solving these equations leads to the following relation describing the onset or nucleation 
condition for superconductivity in the film, 

Considering the extrapolation lengths as material constants imposed by the sample prepa- 
ration (mechanical surface treatment, physical surface deposition technique, or chemical 
modification such as oxidation, ...) the temperature dependence is contained in the variable 
£. In order to obtain direct estimates of the finite-size shift of the film critical point, we 
focus on the following particular cases: similar surfaces and dissimilar surfaces. 

Similar surfaces. 



In this case we assume, for simplicity, b\ = b<i = b and b < 0. We can then work out (IV.3) 
to obtain the critical film thickness L/\b\ as a function of temperature, 

_L = l ln I±iM, (IV .4) 

|6| |6| l-i/W ( ' 

valid for b < and £/|6| < 1 only. This remarkable relation describes the increase of the 
film critical temperature T C (L) upon reduction of the film thickness L. It is graphically 
represented in Fig. 4. Since £/|6| oc \T — T c | _1//2 , and the surface critical point of the semi- 
infinite system (with a single surface) lies at T cs > T c such that £(T cs )/b = — 1, relation 
(IV. 4) predicts that T C (L) > T cs . A convenient temperature variable for our purposes is 
b 2 /£ 2 = (T - T C )/(T CS - T c ), which we denote by t (see Fig.4). With this notation, t c = 
and t^ — 1, as in previous sections. 

Eq.(IV.4) can be written in more compact form if we express the film thickness as L/£, 

|i=tanh| (IV.5) 

but we remark that sample size and temperature are mixed when using the variable L/£. 
We distinguish two regimes: the macroscopic regime L ^> £, and the microscopic regime 
L ^C £. In between is the mesoscopic range L ps ^. In the macroscopic limit (IV.5) gives the 
critical-point shift, 
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£/6 = -l + 2e- L/? , (IV.6) 

or, using the definition t = b 2 /!; 2 , 

t c (L) = 1 + 4e- L/? (IV.7) 

Clearly, for L — ► oo, while keeping the coherence length finite, the ratio £/6 converges 
exponentially fast to — 1. This signifies that T C (L) decreases towards T cs . In contrast, simple 
finite-size scaling ideas would suggest that T C (L) converge to the bulk T c . This is not the 
case for our system, and therefore the finite-size shift in zero field is anomalous. 

For microscopically small thicknesses (L <C £) the increase of the critical temperature 
follows the power-law 

t c (L) »2|6|/L (IV.8) 

Of course, the achievement of very small thicknesses (submicron to nanometer range) is 
subject to practical limitations. 

The increase of T C (L) upon reduction of the film thickness is analogous to an effect that 
has been uncovered in the context of twinning-plane superconductivity [4]. Certain ma- 
terials (Sn and Nb) display enhancement of superconductivity near an internal twinning 
plane. In a situation with closely spaced twinning planes, the enclosed slab of material 
experiences a transition to superconductivity at T C (L) given by the same relation as (IV. 4). 
Twinning-plane superconductivity is a special phenomenon. We would like to emphasize 
that the increase of T C (L) due to confinement is a more general phenomenon and occurs 
for thin films with surface enhancement, regardless of the precise microscopic origin of the 
enhancement. 

Dissimilar surfaces. 



As a first concrete example, we consider enhancement on one surface only and assume that 
the other surface corresponds to a direct contact with vacuum or an insulator. We take 
b\ < and b 2 = oo. In this case the increase of T C (L) still follows a law similar to (IV. 4) but 
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to obtain the same transition temperature increase a further film thickness reduction by a 
factor of 2 is required, 

w\ = w\ hl T^m- < IV ' 9 > 

In a second example we make the second surface unfavourable to superconductivity by 
introducing suppression of the wave function by assuming 62 > and 62 > \bi\- Physically, 
this corresponds to direct contact with a normal metal. The two surfaces are now in com- 
petition and the increase of T C (L) with decreasing film thickness qualitatively still behaves 
as in Fig. 4, but becomes progressively weaker as 62 is decreased. The enhancement effect is 
lost in the antisymmetric limit 62 —*■ —b\. 

4.2. Cylindrical wire. 

For the case of axial symmetry we adopt cylindrical coordinates and write the linearized GL 
equation in zero external field in the familiar Schrodinger equation form 

2m \r dr dr r 2 dcp 2 dz 2 ) 
The coherence length and thus the temperature is related to the "energy" —a through 
H 2 /2m\a\ = £ 2 , with m twice the electron mass. The boundary condition on the cylinder 
surface, at r = R, reads 



■^ (iv.11; 



dip 

dr ., 

r=R 

with b < for surface enhancement of superconductivity. 

Solutions are of the form ip{r,<p,z) = f '(r)e dlt y e lkz . Since we are looking for the lowest 
"energy", we can set k = and / = 0. Indeed, for k 7^ the energy simply increases by 
a positive amount proportional to k 2 . Also, for I 7^ 0, the energy is strictly greater than 
for I = 0. This is due to the strict positivity of the angular momentum contribution, in 
combination with Ritz' theorem, 

_ < ipi I H I ip t > <ipi\H \ipi> <ipi\ L 2 z /2mr 2 | ipi > 

-t^i = 1 = 1 1 1 

<ipi\ipi> <ipi\ipi> <ipi\ipi> 
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< i/ji I H | i/ji > <i) Q \H \j) > _ 
> — j— ; > ; — j— ; = A) (IV.l^J 

<ipi\ipi> <ipo\ipo> 

Here Hq is the first term of the Hamiltonian, associated with the radial kinetic energy. The 
scalar product is defined on the support r G [0, R] for functions that satisfy the bound- 
ary condition (IV. 11). Note that this boundary condition drives the ground state energy 
negative, corresponding to tunneling states with negative kinetic energy. Note that the 
temperature for onset of superconductivity lies above T c , since a > in the ground state. 

After rescaling r by £, it is easily seen that the ground state eigenfunction is the modified 
Bessel function I (r/£), which has the shape of a "hammock" with a smooth minimum on 
the cylinder axis and a maximum on the surface. Application of the boundary condition 
leads to the nucleation condition, 

1 - ^M (IV 13) 

If we make the identification L = 1R we can compare this critical-point shift for the cylin- 
drical wire with that of the thin film. This is done in Fig. 4, with again t = b 2 /^ 2 . Note that 
t is plotted versus 2R/\b\ instead of 2R/C,. With this choice of variables it is understood that 
b is a fixed material constant, so that Fig. 4 presents a temperature versus diameter diagram. 
It is obvious that the increase of T c for cylinders is stronger than for films. This is seen 
dramatically in the experimentally most relevant asymptotic regime of large radius, R ^> £, 
for which we obtain the power law 

4t » 1 - — (IV.14) 

|6| 2R V ; 

This implies a slowly decaying algebraic critical-point shift 

t c (R)Ril + Z/R, (IV.15) 

so that T C (R) — T cs oc 1/R, in contrast with the exponential decay found for the film. 
On the other hand, for small radii, R ^C £, we obtain 

t c (R)^2\b\/R (IV.16) 
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which is similar to (IV. 8). 

4.3. Spherical grain. 

Using similar arguments as for the cylinder one sees that in zero field it suffices to work with 
the radial wave function, which satisfies the differential equation 

- £?£<"> = - aS < IV ' 17) 

with the same boundary condition (IV. 1 1) but now applicable to the surface of the sphere. 
Writing u — rf and using -u(O) = 0, we find u(r) oc sinh(r/£). Application of the boundary 
condition then leads to the nucleation condition 

a\b\ = coth(R/0-£/R (IV.18) 

which results in an increase of T C (R) upon a reduction of R/\b\ shown in Fig. 4. The effect 
is stronger yet for spheres than for cylinders. 

The asymptotic critical-point shift for large R is algebraic, 

£/\b\Kl-£/R (IV.19) 

implying T C (R) — T cs oc 1/R as for cylinders, but with an amplitude larger by a factor of 
2. Note that the mean curvature differs from that of the cylinder by the same factor. The 
effect of confinement is therefore most pronounced for mesoscopic spherical grains. 
Finally, in the microscopic limit R <C £ we find 

t c (R) « 3|6|/i2 (IV.20) 

Recapitulating, we find that for cylinders and spheres the increase of T C (R) is qualita- 
tively stronger than that of T C (L) for films, in view of the 1/R, or "curvature" -dependence 
of the critical-point shift. 

4.4. Comparison between spherical and cubic grains. 
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In going from a planar film to a cylindrical wire and a spherical grain the effective dimen- 
sionality of the system is reduced from 2 to I and 0, respectively. An alternative way of 
reducing the dimensionality is to go from a planar film to a rectangular rod and a cubic 
grain. In this case, however, the surface is not smoothly curved, but displays strong geo- 
metric singularities in the form of sharp edges and corners. In this subsection we calculate 
the critical-point shift for rectangular sample topology. 

Consider a hypercube in n dimensions (n — 1, 2 or 3) of size L n and extend the system 
infinitely in the remaining 3 — n dimensions. The choice n = 1 reproduces the case of the 
planar film, n = 2 the rod with square cross-section and n = 3 the cubic grain. Note that 
the effective dimensionality is d = 3 — n. The GL equation in zero field takes the form 

h 2 

Aip = -aip (IV.21) 

2m 

with a > since T > T c . The boundary conditions on the faces of the hypercube read 

^{x % = 0) 



Xi=0 



■^ X \ =L) (IV.22) 



dip 
dxi 

dip 

9x i Xl =L 

where % — 1, ...,n and b < 0. 

The technique of separation of variables leads in this case to the exact ground state, 
since, as we shall show, the wave function has no nodes. Thus we assume the product form 

iP({ Xi }) = HUM*) (iv.23) 

Solutions to (IV.21) are of the form 

fi( Xi ) = Ae x ^ + Bie-'M* (IV.24) 

The decay lengths & satisfy the constraint 

]T£- 2 = r 2 = 2ma/ft 2 (IV.25) 

Imposing the boundary conditions leads to the equations 
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1+6/W 



1, ...,n (IV.26) 



In order for the wave function to possess no nodes, and therefore to correspond to the exact 
ground state of the Schrodinger problem equivalent to our system of equations, it is necessary 
to consider only the solutions for which £,i/\b\ < 1. Furthermore, for a given ratio L/\b\ this 
solution is non-degenerate and isotropic, so that 

& = SVd, i = l,...,n (IV.27) 

The critical-point shift for L ^> £ is now given by 

t c (L) w n{\ + 2e- L/( ^) (IV.28) 

The exponential decay of the shift T C (L)— T c (oo) turns out to be characteristic for rectangular 
geometry, in contrast with the power-law decay for curved surfaces. Moreover, since the 
condition ^/b = — 1 differs from the condition £/b = — 1, the temperature of onset of 
superconductivity in the thermodynamic limit L — ► oo depends on n in rectangular geometry. 
We find 

t c (oo) = (T c (oo) - T C )/(T CS - T c ) = n, (IV.29) 

with T cs the surface critical temperature for a single planar surface (in a semi-infinite system) . 

It is a paradox that the transition temperature in the thermodynamic limit L — ► oo 
depends on the shape of the sample. However, after some reflection it is clear that in this 
limit the onset of superconductivity is limited to the vicinity of the surface (for n — 1), the 
vicinity of the wedge where two surfaces meet at right angles (for n = 2) and the vicinity 
of the corners (for n — 3), where superconductivity is strongly enhanced. Away from these 
boundaries, in the interior of the system, very little superconductivity should be expected 
just below T c (oo). Therefore, from a physical point of view, the fact that T c (oo) > T c reflects 
boundary superconductivity rather than bulk superconductivity. 

In the microscopic limit (L <^ £) we obtain 
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t c (L) « 2n\b\/L (IV.30) 

which corresponds well to the results for the curved surfaces, if, as usual, we make the 
identification L = 2R. 

V. CONCLUSIONS 

In this paper we have established a close analogy between capillary condensation in fluids 
and the transition from surface superconductivity to mesoscopic sample superconductivity. 
Furthermore, the interplay between capillary condensation, prewetting and wetting, has 
been studied in superconductors which display an interface derealization transition. In the 
limit of strongly type-I superconductors a full analytic description has been given for the 
finite-size effects on the various phase transitions involved. 

We have scrutinized the anomalous critical-point shift in mesoscopic samples in zero 
field, and the standard finite-size scaling for the transition to superconductivity in non-zero 
magnetic field. The critical-point shift in zero field is anomalous in the sense that T C (L) 
or T C (R) converges to T c (oo) > T c instead of T c . Standard finite-size scaling would have 
predicted T c (oo) = T c . However, in a fixed non-zero magnetic field H, no matter how small, 
the transition temperature in the limit L — > oo (or R — ► oo) converges to the temperature 
associated with the bulk critical field H c = H. In other words, in non-zero field the bulk 
two-phase coexistence line is fully restored in the macroscopic or "thermodynamic" limit. 

The anomaly is thus confined to the temperature segment [T c , T cs ] in vanishing field. 
This segment becomes part of the two-phase coexistence line in the limit of large thickness 
or diameter of the sample. As a result, the capillary condensation line displays a corner or 
bend in the vicinity of (T — T c ,H — 0). How this corner develops can be seen in Fig. 5, which 
shows the phase diagram for L/\b\ = 20. Comparing this with Fig.l, which corresponds to 
L/\b\ = 8 we notice that the capillary condensation line approaches the bulk coexistence line 
in non-zero field, but the supercritical temperature segment remains part of the capillary 
condensation line. 
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We have indicated by FSS1 the finite-size shift of the capillary condensation transition 
in the partial wetting regime. Likewise, FSS2 refers to the finite-size shift in the complete 
wetting regime, and FSS3 marks the shift in the supercritical region. These shifts obey the 
analytical laws derived in Section 3. For this large value of L/\b\ the critical-point shift in 
zero field is so (exponentially) small that T C (L) practically coincides with T cs . 

The significance of the zero-field critical temperature in the macroscopic limit T c (oo) 
must be considered with care. For planar, cylindrical and spherical geometry, T c (oo) = T cs . 
This implies that for asymptotically flat surfaces (with vanishing curvature everywhere) su- 
perconductivity starts near the surface. For macroscopic cubes, however, superconductivity 
can start at higher T because it nucleates first in wedges and corners. This is why, for 
square rods and cubes, T c (oo) is progressively increased. This apparent increase of T c in the 
macroscopic limit disappears when the sharp edges are fully rounded. In the mesoscopic and 
microscopic regimes smoothing of corners or other asperities is irrelevant, since shape details 
are then small relative to the scale of £ and cannot be "resolved" . The critical temperature 
is then independent of the details of the shape of the sample. 

From a more practical point of view our most noteworthy result is the significant increase 
of T C (R) for mesoscopic cylindrical or spherical superconductors with surface enhancement. 
The main asymptotic relation between T C (R) and the sample radius can be summarized in 
the form, to leading order in 1/R, 

t^l + 2c£ (V.l) 

where t = (b/£) 2 = (T C (R) —T C )/(T CS —T c ), and c = (I/-R1 + l/i?2)/2 is the mean curvature. 
For cylinders, c = 1/2R. For spheres, c = 1/R. 

We illustrate the theoretically predicted increase of T c by means of a numerical example 
based on the (old) results of Fink and Joiner [3] for a surface-enhanced Ino.9g3Bio.007 foil, 
with T c m 3.5K, and k ~ 0.37. After cold working of the sample surface the critical 
temperature increased to T cs = T c + 0.02K. An order of magnitude estimate suffices, so we 
take a typical value for the coherence length amplitude, £0 ~ IOOOA and obtain £(T CS ) = 
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£o(T cs /T c — l) -1 / 2 ps 1.3yum. Since b = — £(T CS ) we obtain that |6| is of the order of a micron. 

This estimate allows us to predict the further increase of T c due to confinement. For 
example, for a thin film of thickness L = 10/im of the same alloy the predicted T C (L) 
on the basis of (IV.4) is still T c + 0.02K. The confinement effect is imperceptible for this 
thickness. Reducing the thickness further to L — 1/im we obtain T C (L) = T c + 0.05K and 
for L = 0.1/im we find T C (L) = T c + 0.41K. Now, changing the geometry from planar to 
cylindrical or spherical the increase of T c is more striking. For example, for a spherical grain 
of diameter 10/im of the same material, (IV. 15) predicts T C (R) = T C + 0.03K. For 2R = 1/im 
we obtain T C (R) =T C + 0.13K, and for 2R = 0.1 jum we get T C (R) =T C + 1.21K = 4.71K. 

We would like to stress that cold working of the sample surface is only one of the possi- 
ble surface treatments that can lead to surface enhancement of the superconducting order 
parameter. In modern experiments deposition of a thin layer (of thickness less than the 
coherence length) can be performed in a clean, controlled and reproducible way. This layer 
affects the boundary condition, and measurement of the critical temperature in zero field is 
sufficient to establish the sign and, for b < 0, also the magnitude of the surface extrapolation 
length b. 

Further, our results for zero magnetic field are independent of the GL parameter k and 
thus apply to type-I and type-II superconductors alike. In our opinion the results concerning 
the increase of T c are possibly also relevant for surface-enhanced high-T c superconductors. 
The reason for this belief is that in at least two experiments surface enhancement was 
found in high-T c materials. Fang et al. measured an increase of T c of several degrees K 
in macroscopic YBaCuO samples with twinning planes [18], and Schwartzkopf et al. in 
HoBaCuO [19]. Abrikosov and Buzdin [20] invoked the GL theory with the same boundary 
condition, b < 0, to describe this phenomenon. 

In a future publication we intend to report on calculations of global H — T phase dia- 
grams for k > 0. In addition to the features we have discussed in this paper, second-order 
transition lines and tricritical points appear. Also, the film vortex phase becomes stable in 
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a small region of the phase diagram, and gains importance as k is increased. 
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APPENDIX A: CAPILLARY CONDENSATION APPROACHING COMPLETE 

WETTING. 

In this appendix we are concerned with the derivation of Eq.(III.8). The starting point 
is the condition for capillary condensation ^ pw = i cap . This is worked out as follows, 

l P w - leap = ^1sc,n - (H R - 1/2) (L - 2Z)/£ + 5i pw - 5i cap (A.l) 

The term 5i pw represents the free energy of two constrained wetting layers minus that 
of two equilibrium wetting layers. If we adopt the notation 

J(V>; A, B) = A ' r/2 (A.2) 

l-^ 2 + ij 4 /2 + B 



we obtain 

5l pw = 2 dipJ(ip',12,Hl)-2 # J V; 1/2, 1/2) (A.3) 

JO JO 

The fact that the argument A of the function J equals 1/2 in both terms signifies that Hr is 
set equal to H Rc . The constraint is imposed by setting the argument B equal to H R > 1/2 
in the first term only, so that the constrained layer has the same thickness and profile as 
an equilibrium layer in a field H R . For more details on the physics of constrained surface 
sheaths, we refer to [9]. 

The term 5i cap represents the free energy of the capillary condensed state minus 27^5^, 



31 



hca P = 2 #J(^;l/2,C)-2/ #J(^;l/2,l/2), (A.4) 

where C = C(L) is the "constant" in the first integral for superconducting film states, (II. 3). 
We now examine the calculation of these terms, and begin with S , f pw . We define e = 
H\ — 1/2 and consider an expansion to first order in e. The first contribution comes from 
the e-dependent upper limit of the integral. Explicitly, 

iMO) 2 = 1 + {i/bf + v /(l + (e/6) 2 ) 2 -(l + 2e) (A.5) 

For e — > 0, VV"(0) reduces to ip(0). The second contribution comes from the e-dependence 
of the integrand J(ip; 1/2, H^). After some elementary algebra, we find to first order, 

**" " m2 + l + lim ( f m - f 1 ) # 7 1+ f (A.6) 



v^e 2^(0)(^(0) 2 -l-(e/6) 2 ) ^Vi Vo ; r (l-^ 2 ) 2 + e 

The evaluation of the limit must be done carefully. We have separated the integrand 
into a symmetric and an antisymmetric part (with respect to ip = 1). The symmetric part, 
which contains the dominant singularity proportional to l/(ip — l) 2 , is largely canceled by 
subtracting the integrals. In the part that remains one may safely set e = and the result 
is the contribution, with x = ip — 1, 

f 1 8 - 2x 2 + x 4 

6sym ~ " i(o)-i dX * 2 (4-* 2 ) 2 (AJ) 

which is simple to evaluate. In the antisymmetric part, however, one may not exchange the 
limit e — *■ with the integration. One must calculate 

/ t-\ rip(o)-i\ 2ex — 2x 5 

S an u sym = Jim (^ +J Q jdx e2 + 8ex2 + 2ex4 + 1Qx4 _ 8x , + x s (A.8) 

This leads to two contributions, associated with the two terms in the numerator of the 
integrand. The result is, with z = x 2 /e, 



2 ( r 1 /-V>(0)-i\ 2x 

o "~{i + AzY ~ [Jo + h ) dx j4^x~ 2 y 



SanUsym = I dz - ( I + dx (A.9) 



which is elementary to evaluate. 
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The final result is 

8'y pw = y/2e + o{€) (A. 10) 

The correction o(e) goes to zero faster than e, in the manner ae 2 log(l/e), according to 
numerical computations. For example, for £/6 = —1 we obtain a m 0.26. It is interesting 
to note that the leading-order term is independent of ip(0), and thus independent of the 
"temperature" variable £/&. 

We now turn to the evaluation of <57 cap . As a first step we examine the dependence 
C(L) for capillary condensation states more closely. For this it suffices to study the leading 
terms in L/C, for C — ► 1/2 from below. Using (II. 4), together with the relations ifj m = 
(1 + (1 - 2C) l l 2 ) 1 ' 2 and ^ cap (0) 2 = 1 + (£/b) 2 + ((1 + {C/bff - 2C) 1/2 , we obtain after some 
algebra, 

^7f ln I75b + D + o(1) (A ' n) 

where o(l) denotes terms that vanish for C — ► 1/2. We verified numerically that the 
constant £) is typically of order unity and depends on £/6. This result implies the (expected) 
exponential decay of C{L) for large I, 

(1/2 -C) oc e -^ L /« (A.12) 

In a second step we find, using numerical computation, that for C approaching 1/2, 

hca^-^l-C) (A.13) 

so that the finite-size correction for the surface free energy of the capillary condensed profile 
decays exponentially rapidly with the thickness of the film, in view of (A.12). The fact that 
^7cap is negative expresses the lowering of the free energy of the superconducting film by 
confinement. This is a manifestation of the fairly general observation that "like" surfaces 
attract each other [21]. In contrast, ^7 pw > 0, reflecting the repulsion between the SC/N 
interface and the surface. 
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FIGURE CAPTIONS 

Figure 1. 

Capillary condensation phase diagram for strongly type-I superconductors with surface en- 
hancement, for a film of thickness L/\b\— 8, in units of the surface extrapolation length b 
(b < 0). The magnetic field H is scaled with the field Hp at the interface derealization 
transition D. The temperature variable is t = (T — T C )/(T CS — T c ). The capillary conden- 
sation transition runs mostly parallel to the bulk coexistence line. It meets the prewetting 
transition at a film triple point Tl, and again at T2. It terminates at the film critical point 
in zero field, at T C (L) > T cs , but very close to T cs (imperceptible difference). The solid 
lines indicate first-order phase transitions for the film. The dashed line is the (metastable) 
continuation of the prewetting line. 

Figure 2. 

Three coexisting film phases at triple point Tl of Fig.l. A normal phase (ip = 0) coexists 
with a surface superconducting film (with two sheaths), and with a capillary condensed 
superconducting film. 

Figure 3. 

Three coexisting film phases at triple point T2 of Fig.l. Under these conditions the sample 
in bulk would show no superconductivity at all. 

Figure 4. 

Increase of the critical temperature as a function of the thickness or diameter of meso- 
scopic superconductors with surface enhancement. The temperature variable is t c (L) = 
(T C (L) — T c )/(T c (oo) — T c ), where L is the film thickness, to be replaced by 2R for cylinders 
or spheres. T c (oo) equals the surface critical temperature T cs of a semi-infinite sample, and 
exceeds T c . The important difference between curved and planar surfaces is the long-range 
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(algebraic) decay of the critical-point shift for large radius. 

Figure 5. 

Finite-size scaling (FSS) of the capillary condensation transition at k — illustrated for 
a thick film, with L/\b\ = 20. See also Fig.l for comparison. The arrows FSS1 show the 
algebraic shift of order 1/L between the film transition and the bulk transition, in the par- 
tial wetting regime (no surface sheaths intervene). FSS2 indicates the similar shift in the 
complete wetting regime (with surface sheaths induced by the prewetting transition). For 
T > T c the exponentially small shift is apparent (FSS3) and clarifies how the capillary con- 
densation line eventually converges, for large L, to a corner shape or "dog leg" consisting of 
the bulk coexistence line supplemented with the segment [T c , T cs ] on the temperature axis. 



37 



